Pemrograman Prosesor Paralel Besar: Pendekatan Praktis: Pengantar Visualisasi Molekuler dan Peta Potensial Elektrostatik

Visualisasi molekuler berfungsi sebagai jembatan penting antara koordinat atom dan intuisi biologis. Dengan memanfaatkan perangkat lunak seperti Visual Molecular Dynamics (VMD), para peneliti dapat mengubah data numerik mentah menjadi lingkungan 3D interaktif yang mengungkap tata gerak struktural kehidupan.

1. Peta Potensial Elektrostatik

Peta Potensial Elektrostatik adalah representasi berbasis kisi 3D yang menunjukkan distribusi muatan listrik di seluruh molekul. Setiap voxel dalam kisi menghitung jumlah potensi listrik dari semua atom: $$V_j = \sum_{i} \frac{q_i}{r_{ij}}$$. Peta-peta ini berfungsi sebagai pengganti medan gaya, mengidentifikasi daerah dengan afinitas tinggi untuk pengikatan dan lipatan.

2. Keunggulan GPU

Menghitung peta-peta ini sangat mahal secara komputasi. Seperti yang ditunjukkan pada Gambar 9.1, proses tersebut melibatkan rendering pita protein kompleks yang dikelilingi oleh awan titik padat berwarna (merah untuk negatif, biru untuk positif). Paralelisme besar ini membuat GPU sangat ideal untuk simulasi semacam ini.

3. Penjumlahan Coulomb Langsung (DCS)

DCS adalah algoritma pilihan untuk pembuatan peta. Algoritma ini bergantung pada instruksi rsqrtf untuk perhitungan akar kuadrat balikan berkinerja tinggi, memanfaatkan memori konstan untuk menyebarkan data atom ke semua thread pemrosesan secara bersamaan.

TERMINALbash — 80x24

> Ready. Click "Run" to execute.

QUESTION 1

Compare the number of operations (memory loads, floating-point arithmetic, branches) executed in each iteration of the kernel shown in Figure 9.7 compared to that in Figure 9.5. Keep in mind that each iteration of the former corresponds to four iterations of the latter.

Figure 9.7 increases memory loads by 400% compared to 9.5.

Figure 9.7 reduces branch instructions by 75% and reuses atom coordinate loads for 4 grid points.

Figure 9.5 is more efficient because it uses fewer registers per thread.

There is no difference in the number of operations when normalized by grid point.

QUESTION 2

What are two potential disadvantages associated with increasing the amount of work done in each CUDA thread, as shown in Section 9.3?

Lower global memory bandwidth and increased warp divergence.

Increased register pressure and potential for reduced hardware utilization/occupancy.

Reduced CPU-GPU transfer speeds and higher thermal throttling.

Increased constant memory size requirements and cache misses.

QUESTION 3

Why is 'Constant Memory' utilized for atom data in the DCS kernel?

Because it is faster than registers for local storage.

Because all threads in a warp access the same atom data at the same time, benefiting from a broadcast.

Because it allows for atomic write operations.

Because it bypasses the L1 cache entirely.

QUESTION 4

In the context of VMD, what does a red area on an electrostatic map typically represent?

A region of high positive potential.

A region of high negative potential.

A hydrophobic lipid boundary.

An area where the GPU has failed to calculate data.

QUESTION 5

Which CUDA function is prioritized for distance-based potential calculations in DCS?

sqrtf()

rsqrtf()

powf(x, -0.5)

fast_div()